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Abstract. This paper presents a numerical investigation of isolated filament 
dynamics in a simulation geometry representative of the scrape-off layer (SOL) of 
the Mega Amp Spherical Tokamak (MAST) previously studied in [N.R.Walkden et.al^ 
Plasma Phys. Control. Fusion, 55 (2013) 105005]. This paper focuses on the evolution 
of filament cross-sections at the outboard midplane and investigates the scaling of the 
centre of mass velocity of the filament cross-section with filament width and electron 
temperature. 

By decoupling the vorticity equation into even and odd parity components about the 
centre of the filament in the bi-normal direction parallel density gradients are shown to 
drive large velocities in the bi-normal (approximately poloidal) direction which scale 
linearly with electron temperature. In this respect increasing the electron temperature 
causes a departure of the filament dynamics from 2D behaviours. 

Despite the strong impact of 3D effects the radial motion of the filament is shown to be 
relatively well predicted by 2D scalings. The radial velocity is found to scale positively 
with both electron temperature and cross-sectional width, suggesting an inertially 
limited nature. Comparison with the two-region model [J. R. Myra et.al^ Phys. 
Plasmas, 13 (2006) 112502] achieves reasonable agreement when using a corrected 
parallel connection length due to the neglect of diamagnetic currents driven in the 
divertor region of the filament. 

Analysis of the transport of particles due to the motion of the filament shows that the 
background temperature has a weak overall effect on the radial particle fiux whilst the 
filament width has a strong effect. 
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Controlling particle and heat fluxes to material surfaces is an essential step towards 
achieving net energy gain from a high powered fusion device. Excess heat fluxes to 
divertor surfaces will result in melting damage requiring repair and limiting operation. 
Excess particle fluxes can cause erosion of materials, again causing damage but also pro¬ 
viding an influx of impurities into the plasma thereby degrading its performance. In the 
last decade it has become clear that particle transport into the scrape-off layer (SOL) 
cannot be accurately captured by local transport models such as the diffusion-advection 
paradigm [T] . Instead over 50% of particle transport can be mediated by the propagation 
of large, coherent, fleld aligned structures called filaments (also commonly referred to as 
blobs) [2]. Filaments have been observed on many tokamaks [SI 51 El El 0 El El Uni EH 
and other magnetically confined plasma devices da EH ES] suggesting a rather uni¬ 
versal nature. Recently the formation of a density shoulder in the SOL density profile 
dH EH EH ED has been linked to an increase in filamentary transport jlS] the gas 
fuelling rate of the plasma is increased. This may be favourable as it spreads the impact 
of erosion over a larger area on the divertor surface. It runs the risk however of causing 
excess interaction with first-wall materials, leading to damage of the first wall and dilu¬ 
tion of the core due to an influx of first-wall impurities [20] . Predictive control of these 
effects is desirable, however before any level of prediction can be made, a fundamental 
understanding of the processes governing filament motion must first be established. 

Models of filament motion have been traditionally based on 2D interchange motion 
[m EH EH El EH [26] where variation in the magnetic field structure leads to a non-zero 
divergence of the diamagnetic drift of ions and electrons. This arises due to the inho¬ 
mogeneity of the background magnetic fleld and drives a diamagnetic current across the 
filament. This drive is balanced by dissipation either due to aerodynamics drag, termed 
the inertially limited regime, or by sheath currents, termed the sheath limited regime. 
In each case the result is that the filament propagates radially due to an E x B velocity 
as a result of a dipole potential forming across its cross-section. The formation of the 
potential dipole has been observed experimentally [2ZIEHI- For an excellent review of 
2D models of filament motion and their comparison with experiment see ref [221 • In 
recent years models for filament motion have been extended to 3D in a basic slab geom- 
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etry [301 [an El and in more realistic geometries [aaiiaiiias]. This paper follows from 
this work and investigates in detail the motion of filaments at the outboard midplane 
and how changes in temperature and filament width affect the resulting transport of 
particles. 

This paper is organized as follows: Section introduces the model used for filament 
simulations and describes the role played by 3D effects. Section describes briefly the 
magnetic geometry employed to simulate the MAST SOL, as well as providing details 
of the numerical methods used in BOUT++ [3S]. Section presents the dynamics of 
filaments across a temperature and filament width scan at the outboard midplane and 
compares the results to the two-region model [HU EH]. Section gives a detailed analysis 
of particle transport due to filament motion and assesses dominant factors. Finally 
discusses the validity and implications of the results and concludes. 

2. Model 

The model used in this paper to investigate the dynamics of filaments is based on 
the continuity of density and current. The model follows the drift ordering [HH] and 
additionally makes the Boussinesq approximation and assumes that the plasma is 
isothermal and electrostatic. This prevents the model from capturing effects due to 
magnetic or thermal fluctuations. The low pressures in L-mode filaments lead to low 
values of /5, however high pressure gradients may still make the electrostatic assumption 
fragile [iniiii]. Since this work seeks to establish the connection between basic models of 
filament motion and 3D models these effects have been neglected, however their inclusion 
should be treated as a priority for future study. Since thermal conduction time scales 
are much faster than particle loss time scales (typically L||/cs), temperature fluctuations 
are likely to be drained to the divertor much more quickly than density fluctuations, so 
the isothermal assumption may be justified in the far SOL. This is a fragile assumption 
however and non-isothermal effects will be considered in a future paper. For the present 
paper it is instructive to use a highly reduced a model to isolate the causes of any effects 
in the filament dynamics observed. In this vein ions are assumed cold such that T* Tg 
which eliminates FLR effects associated with the diamagnetic contribution to vorticity 
|121I1H|. Finally the effects of parallel ion streaming are neglected as slow in comparison 
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to the cross-field evolution of the filament and electron inertia is neglected. These latter 


two assumptions have recently been shown to be justified in 3D slab simulations [32] • 

Under this set of assumptions the governing equations of the system are: 



(jri 1 

—+ V-(nv) = 0 ^ — = 2c,p4-Vn+-V||J|| 

(1) 


V-J = 0 ^ dt = g^ll'^ll 

(2) 

with 

V\\J\\= TeVii (ln(n) - 0) 

(3) 


d d ^ 

(4) 

and 

ve = CsPsh X 

(5) 


where n is the electron density (ion density is inferred through quasineutrality), J 
is the total current density, J|| is the component of current density parallel to the 
magnetic field and 0 = is the potential normalized to the electron temperature 

Tg in eV. ^ = (V x (b/B))B ^ b X Av is the polarization vector arrising from 
the divergence of the magnetization drift where b is the magnetic field tangency 
vector and k. = b ■ Vb = — b x V x b is the magnetic field curvature vector. 
Cg = ^JeTf,/mi is the Bohm sound speed, ps = Cg/kli is the Bohm gyro-radius and 
fli = eB/mj is the ion gyro-frequency, r^n is the parallel resistivity here taken as the 
Spitzer resistivity, r]\\ = mg/(1.96ne^rg) with the electron ion collision time given by 
Tg = Sy/m^Tg/^/(4-\/^nAe^) = 3.44 x 10^^ (e^) ^ where A ~ 10 is the 

Coulomb logarithm. Terms on the RHS of equation [T] are smaller than leading order, 
however Angus et.al have shown that they must be retained to properly capture the k± 
dependence in the linear regime of the resistive drift wave instability [301 EH- 
Three-dimensional effects are incorporated into the system through parallel Ohm’s Law, 
equation The inclusion of these effects introduces a non-dipolar component of the 
electrostatic potential, thereby inducing motion in the poloidal direction alongside the 
filament’s radial motion [HU [331 132] • To illustrate this the symmetry properties of 
the filament potential can be considered. In the idealized situation, often studied in 
two and three dimensions, the density is taken as a monopolar function with reflective 
symmetry in both the normal (radial at the midplane) and bi-normal (approximately 
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poloidal) directions with the origin taken at the centre of the filament, where the density 
is initialized as a Gaussian function. The density and potential can be split into 
components with an even and odd parity about the filament centre in the bi-normal 
direction. These are labelled as and for the even components and n~ and (/)“ for 
the odd components. The density is assumed to take the form of an even function such 
that n = with n~ = 0. The electrostatic potential is a superposition of both parity 
components, ({) = (j)'^ -\- (j)~. This deconstruction is described schematically in figure [1} 
Substituting the expressions for n and in terms of n+, (/)+ and (/)“, into equation and 


Even Odd 

component component 


Total 



Figure 1. Schematic illustration of the deconstruction of the potential and density 
profiles into odd and even functions about the filament centre in the bi-normal 
direction. The density is assumed to be entirely even whilst the potential is a 
superposition of the two components. 


noting that even order differential operators in the poloidal direction retain the parity 
of the function being operated on whilst odd order operators reverse the parity, gives 

FF In (n) - = (6) 

+ C 3 p 4 -Vln (n) - 

n " 

where the odd and even parity components of the E x B velocity are 

= CsPsh X V 0 + ( 7 ) 

Vg = CsPsh X 

In equation all even parity terms have been placed on the LHS whilst all odd 
parity terms are on the RHS. For equation to be generally true both sides must be 
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independently equal to 0. This allows equation to be split into two coupled equations 
for the odd and even components of the potential. These are 


pI 


dV 


2 A- 


ve ■ vvi(/)+ + v+ • vvi4> 


2 A- 


dt 


= C3p4-Vln(n)-^V;jr(8) 


n 


for 0“ and 

= ^ In (n) — (9) 

for The inclusion of the V||n term in parallel Ohm’s law provides a source for 
the even component of potential. It is commonplace in 2D theories of blob motion to 
neglect this term [211 EHl 1221 HU ESI EH] , though it is possible to accomodate it within 
2D models [l5l 06]. If the electrostatic potential is initialized as 0 or an odd function, 
again commonplace in 2D simulations, then equation prevents any growth of the even 
component in the absence of the V||n term and the potential can only adopt an odd 
symmetry; in this sense simulations neglecting V||n are a rather specihc class of models 
which do not permit a breaking of the odd symmetry of the potential (without the 
addition of extra physics beyond the basic models). It should be noted that if the 
even component of the potential is of a similar magnitude to the odd component then 
significant motion of the filament can occur away from the radial direction. This violates 
the assumption that n~ = 0 and the present analysis breaks down. Never the less this 
serves to illustrate in a general manner the role played by 3D effects in the motion of 
filaments. By balancing the source for the odd component against the source for the even 
component a critical temperature can be approximated above which the initial growth 
of the even component is comparable to the odd component and 3D effects may impact 
the motion of the filament. Taking V|| ~ l/-b|| and V_l ~ 1/(5x as approximations where 
L|| is the parallel length scale of the filament and is the perpendicular radius of the 
filament, the critical electron temperature is estimated as 


pI 


dVlq 

dt 


+ Vr-VV^(/> +vt-VVi 


2csPs\$\ln{n) Tecr|| 




nU 


ln(n) ^ Te 


2nmeL?, |^| 


2/3 


1.96 X 3 X W05±eB 


4el^(10) 


where |^| = 1/R with the major radius R = 1.5m, = 2cm, n = 5 x and 

B = 0.3T are taken as representative values of the MAST SOL [1^, whilst L\\ = 8m is 
taken directly from the simulation geometry used here (see section . The numerical 


coefficients in the denominator of equation 10 arise from the definition of and rj\\. 


This serves as an approximate temperature above which 3D effects can be expected to 
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play a role in filament dynamics. Given that near SOL temperatures (and in some cases 
far SOL temperatures) in MAST are universally above 4eV |lZj it is important to assess 
the role of these 3D effects. To investigate the role of 3D effects on filament motion and 
the resulting particle transport scans in both electron temperature and filament radius 
have been conducted in a geometry representative of the SOL in MAST. 

3. MAST Geometry 

The simulations in this paper have been conducted within the BOUT++ framework 
[36l 08]. The simulation geometry is a flux-tube in the SOL of MAST constructed 
around a field line on the = 1.15 flux surface where — 'ipsep) / {'4’axis — i^sep) 

is the normalized poloidal magnetic flux and 'tp, ipaxis and prsep are the poloidal magnetic 
flux on the flux surface, magnetic axis and separatrix respectively, ippp = 1.15 gives a flux 
surface suitably far from the separatrix to consider filaments as isolated structures, but 
sufficiently far from the first wall to study their motion. Figure shows the ippp = 1.15 
flux surface in MAST, as well as a 3D projection of the flux-tube domain. Limitations 


2.0 1.5 1.0 0.5 0.5 1.0 1.5 2.0 

R (m) R (m) 


Figure 2. Left: ipN = 1-15 flux surface in a MAST connected double-null Ohmic 
L-mode plasma. Right: 3D projection of the flux-tube domain constructed around a 
field line on the ipN = 1-15 surface 

of the fiux-tube approach used here are discussed in [HS]- 

The simulations are carried out in the field-aligned ballooning coordinate system defined 
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by the coordinates 


x = 'ip-'ipo y = e z = C - dO' 


( 11 ) 


Fourth order central differencing is applied to diffusive operators in the x and y 
dimensions whilst a 3rd order WENO (Weighted Essentially Non-Oscillatory) scheme is 
used for parallel advective operators [19]. EETs are used in the angular coordinate z. 
Neumann boundary conditions are applied to all variables on the boundaries in x whilst 
periodic boundary conditions are applied in z with a toroidal period of 15. This allows 
filaments to be simulated with a toroidal mode number of 15, which is on the lower 
bound of the experimentally observed filament quasi-periodicity nDj. This periodicity 
was chosen as it allowed filaments to be treated as isolated; their profiles in density or 
potential did not reach the periodic boundary, particularly as the potential dipole was 
established. The y coordinate labels the position along the magnetic field line. Sheath 
boundary conditions are applied to the boundaries in y by setting the parallel current 
J\\ to the current at the entrance to the magnetic pre-sheath Jgh isniEi] where 


J% = ±.necs (l — exp (pw — ') 


( 12 ) 


with the wall potential, pw, set to 0 since it acts as a constant reference in the isother¬ 
mal case. See |36l |33j for more details concerning the coordinate system. The system is 
integrating in time using the PVODE [52] implicit time-solver which is a Jacobian-free 
Newton-Krylov solver using the GMRES Krylov subspace method. 512 grid points have 
been used in x, 256 are used in z and 64 are used in the parallel y coordinate. This 
provides sufficient resolution to capture the most unstable resistive drift-wave when dis¬ 
sipation is removed from the system [35]. Eilaments are initialized with a 2D Gaussian 
density in the perpendicular plane with a peak density of twice the background. They 
are initially homogenous along the magnetic field line. The potential is initialized to 
zero and grows in the initial phase of the simulation. Unlike ref [32] the background is 
flat. 


As illustrated in ref. [33] the magnetic geometry of the SOL gives rise intrinsically 
to gradients in the filament density along the magnetic field. There are two primary 
causes of this: 1. Ballooning of the filament near the midplane [SH] and 2. Enhanced 
dissipation in the divertor region [35] . Both of these effects are a result of the variation in 
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the magnetic geometry of the flux-tube along the field line. In the divertor the shearing 
and contraction of the flux tube in the bi-normal direction lead to enhanced cross-field 
density gradients. As a result the diffusion of density is strongly enhanced in this region, 
even at dissipation levels well below that which may be expected experimentally [53] . 
By providing mechanisms for the formation of parallel density gradients the magnetic 
geometry naturally drives the source for the even component of the potential. 

4. Dynamics of the filament cross-section 

4 ..I. Effect of the electron temperature 

Figure]^ shows the evolution of the filament cross-section at the midplane with constant 
background temperatures Tg = leV and Tg = 20eV. The dynamics of the filament vary 




Figure 3. Top: Filament cross-section in the perpendicular plane at the outboard 
midplane with Tg = leV. Bottom: As top but with Tg = 20eV. Colour represents 
density and contours represent electrostatic potential. Increasing the electron 
temperature increases the prominence of the Boltzmann source which leads to a 
departure from interchange dynamics. 
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dramatically between the two cases. At leV the behaviour of the filament is reminiscent 
of the inertially limited regime found in 2D models where a dominantly dipolar po¬ 
tential causes the filament to ’mushroom’ outwards. Resistive drift-waves are observed 
to form on the filament front during the latter part of the filament’s motion at leV, 
as described by Angus et.al [301 EH- These instabilities, shown in figure form when 
density gradients on the filament front steepen as a result of its motion, and have a 
typical length scale much smaller than that of the filament cross-section. 

0.12 
0.10 
0.08 
^ 0.06 
0.04 

0.02 
0.00 

0.00 0.02 0.04 0.06 0.08 0.10 0.12 
x(m) 

Figure 4. Filament cross-section at the outboard midplane with Tg = leV 113/is after 
initialization. Resistive drift-waves form on the front of the filament due to density 
gradient build up, but do not significantly affect the filament. Their length scale is 
much smaller than that of the filament. 



At 20eV the dynamics of the filament depart from the leV case. The dipole poten¬ 
tial is not sustained, instead showing an initial rotation before a break down towards 
smaller length scales. The rotation of the dipole leads to motion in the poloidal direc¬ 
tion, similar to that observed in refs. jsniis], before the break down of the potential 
acts to disperse the filament. To show that these features are a result of the inclusion 
of the V|In term in equation the 20eV simulation has been re-run with this term now 
excluded. Figure [^compares the filament after 12/is of evolution with and without the 
V||n term. Removing the V||n term allows the filament to return to the dipole potential 
by removing any source for an even component of the potential. 

Increasing temperature decreases the resistivity of the plasma by (thus a 20 times 

increase in Tg corresponds to a 90 times decrease in resistivity). Equation suggests 
that as resistivity is decreased the filament should adopt a Boltzmann response such 
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Figure 5. Top; Filament cross-section at the outboard midplane with (left) and 
without (right) the Boltzmann source included. Bottom; V|| ln(n) (colour) and 
V||^ (contours) for the cases with and without the Boltzmann source included. The 
Boltzmann source drives the filament towards a Boltzmann response and is responsible 
for the observed dipole rotation. 


that V||ln(n) ~ V||0. This has been tested in the bottom frames of figure and, as 
expected, the case where V||n is included in the system does indeed tend towards a 
Boltzmann response, whilst the counter case does not. Notably this shows that 2D 
models which inherently neglect the V||n term are not compatible with the requirement 
for the plasma to obtain a Boltzmann response, even at relatively modest temperatures 
such as 20eV. As a result care should be taken when applying 2D models to situations 
where even modest parallel density gradients can occur. 


The even and odd parity components of the potential, 0+ and (f) can be extracted 
from the full potential solution, 0, using 

1 (</■(- - ^o) + '#' (zci - z)) 


(13) 
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0 =\{(I>{z-zq)-^{zo-z)) (14) 

where z is the bi-normal coordinate and Zq is the centre of the filament in the perpen¬ 
dicular plane. Figure shows the evolution of the absolute maximum of the odd and 
even components of the potential during the simulations of the leV and 20eV filaments. 



Figure 6. Absolute maximum of the even and odd components of potential 
(normalized to the absolute maximum of the total potential over the entire simulation) 
during the evolution of the leV and the 20eV filaments. In the leV case the growth 
of the odd component is significantly faster than the even. In the 20eV case the two 
are comparable. 


The relative growth of the even component of the potential increases dramatically with 
increasing electron temperature due to the increased strength of the Boltzmann source. 
In the 20eV case the poloidal motion breaks the even symmetry of the density early 
in the simulation, which causes the arguments based on parity to break down. The 
initial growth of the two components of potential occurs before this symmetry breaking 
though and confirms that increasing electron temperature increases the drive for the 
even component of potential. 


A useful quantity for assessment of the filament motion is the velocity in the nor¬ 
mal {Vip/ IV'^I) and bi-normal (b x Vip/ \h x V'0|) directions respectively. These are 
calculated by projecting the E x B velocity onto the normal and bi-normal unit vectors 
to give 




— Ve ' • 


RBq 


(15) 


and 
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VE-eg = VE- {'V'lp xh) =ve 


RBf) 


(Vz + IF/ij) 
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(16) 
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where is the velocity in the normal direction and is the velocity in the bi¬ 
normal direction, R, Bg and I are the major radius, poloidal magnetic field strength 
and integrated magnetic shear of the flux tube. The convention adopted in ref [Slj 
where the bi-normal direction is labelled with g is adopted here. The centre of mass 
velocity of the filament can then be calculated from 


^e’^) ~ J J ^E’^SndxNdXg/ J Sndxjsfdx^ 


(17) 


where 5n is the density of the filament on top of the background. Figure shows the 
evolution of (v^) and (v^) for each filament in the temperature scan with Sn/no ~ 2 
and fix = 4cm. The velocity in each case shows an initial acceleration period as the 



Figure 7. Normal (left) and bi-normal (right) centre of mass velocity of the filament 
cross-section at the outboard midplane over the scan in Tg. 


potential grows, followed by a peak and a deceleration period. In both the normal 
and bi-normal direction the maximum speed attained by the filament increases with 
electron temperature. In figure]^ the maximum speed is given as a function of electron 
temperature. Figure [^indicates that (v^) oc whilst {\v^\) oc Tg. The scaling of the 
normal velocity is consistent with the inertial regime of filament motion [26] where the 
characteristic velocity has the scaling 

- c,^2|^|fix oc (18) 

A scaling for the bi-normal velocity can be derived by balancing order of magnitude 
terms in equation]^ If the vorticity term, which acts to mix even and odd components 
of the potential, is balanced against the drive term and the odd component of the 
potential is taken to be inertially limited (as suggested by the present simulations) the 
binormal velocity obtains the scaling 

Cs(T\\B 


5 
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Figure 8. Maximal velocity in the normal (left) and bi-normal (right) direction. Inset 

1/2 

is the maximal normal velocity as a function of Te ^ demonstrating a reasonably linear 
relationship. 

On the other hand if the drive term is ballanced by the parallel dissipation term such 
that the Boltzmann response is satisfied the bi-normal velocity obtains the scaling 

^ In ((5n) oc (20) 

The trend of the maximal bi-normal velocity observed over the temperature scan, shown 
in hgure [8| indicates that the Boltzmann response scaling is appropriate to describe the 
filaments observed here. 


Despite the increase in the peak velocity, a decrease in the time taken to reach the 
velocity peak and a speed-up in the deceleration following the peak is observed. These 
factors all affect the net displacement of the filament and may cause it to scale differ¬ 
ently to the velocity. Figure]^ shows the net displacement of the filament centre of mass 
in the normal and bi-normal directions. 


The linear increase in the bi-normal velocity with temperature translates into a linear 
increase in the bi-normal displacement of the filament by ~ 5cm across the temper¬ 
ature scan. By contrast the increase in the normal displacement is weaker than the 
corresponding increase in velocity. This is because, although the peak velocity of the 
filament increases with temperature, the deceleration after the peak is also quicker as 
the break down into drift-wave turbulence occurs. As a result a much more modest 
~ 0.6cm increase in the normal displacement occurs across the temperature scan. 
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Figure 9. Net displacement of the filament cross-section at the outboard midplane in 
the normal (left) and bi-normal (right) direction calculated by integrating the velocity 
curves in figure over the simulation time. 


4-2. Filament Width Scan 


To compliment the scan in electron temperature and further investigate the dynamics 
of filaments in MAST a scan in filament width has been conducted in the range 2cm 
to 8cm at a constant electron temperature of Tg = 5eV , which is characteristic of the 
far SOL of MAST [55l HT], and with Sn/no ^ 2 held constant. This encompasses the 
range of filament widths found experimentally on MAST pUl ITT] . The quoted filament 
widths here refer to the width in the bi-normal direction at initialisation where the 
density exceeds a 1% increase over the background (this is also the limit of the colour 


levels shown in any cross-section image of the filament). Figure 10 shows the midplane 
cross-section of filaments in the S± scan. 
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Figure 10. Filament cross-sections sampled 38/is after initialization for varying values 
of the initial cross-sectional width (in the bi-normal direction) The prevelance of 
3D effects reduces as the filament width increases from left to right. 


The impact of 3D effects is stronger in smaller filaments, where the density shows 
strong bi-normal motion and the potential begins to break down into smaller scale 
structures. By contrast the larger filaments show a dominantly dipolar potential which 
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leads predominantly radial motion. This suggests that, for the parameters used here, 
larger filaments move mainly in the radial direction, whilst smaller filaments can have a 
significant bi-normal component to their motion. Figure shows the evolution of the 
even and odd components of the potential at the limits of the scan. 


— 2cm <5_|_ = 8cm 



Figure 11. Growth of the absolute maximum of the even and odd components of 
the potential at = 2cm and = 8cm. The growth of the two components is 
comparable in the smaller filament but the growth of the odd component is faster in 
the larger filament. 


The relative growth of the even component of the potential decreases as the filament 
width increases. This occurs because at larger filament widths dissipation due to in¬ 
ertia diminishes quicker than the drive due to the diamagnetic current reduces. As a 
result there is comparatively more current flowing across the filament in larger filaments 
which must be closed by current flow along the magnetic field line, thereby increasing 
J||. Since the resistivity, 77 , is constant Ohm’s law (equation implies that an increase 
in the parallel current will drive the filament away from the Boltzmann response. This 
then leads to the observed decrease in the relative strength of with respect to 0 “ as 
fix is increased. 

Figure [T^ shows the centre of mass velocity of the filament at the midplane in the 
normal and bi-normal directions in a subset of simulations contained in the scan. 

The maximum velocity of the filament decreases weakly with increasing in the normal 
direction, in contrast to the scaling predicted in the inertially limited regime. The 
bi-normal velocity also shows significantly more complex behaviour than the scan. 
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Figure 12. Centre of mass velocity in the normal (left) and bi-normal (right) 
directions over the scan in 6±. 


Whilst the peak velocity does decrease with increasing as predicted by the Boltz¬ 
mann response scaling, it does not show a clear dependence. This is likely due to 
the fact that the mechanism for bi-normal motion is more complex than for the normal 
motion and involves the break down of the potential into small scale structures. As a 
result this may not be particularly well represented by a characteristic velocity based 
on an order of magnitude scaling of a ballistic velocity. In fact the even component 
of potential cannot, on its own, lead to non-circnlatory motion and requires some non¬ 
linear interaction with the odd component to produce bi-normal motion. Such complex 
motion is unlikely to adhere to simple scaling arguments. 


The net displacement of the filament is presented in figure 13 


The displacement in 




Figure 13. Net displacement of the filament in the normal (left) and bi-normal (right) 
direction across the 6± scan measured by integrating the curves shown in figure over 
the simulation time. Inset is the net displacement in the normal direction as a function 
of 6Y‘^ displaying an approximately linear relationship. 
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1 /2 

the normal direction obtains a scaling of 5^ suggesting that the inertially limited 


regime may be well represented by the average velocity of the filament over its lifetime 
(i.e. displacement / simulation time). In the bi-normal direction the displacement, as 
with the velocity, shows a general decrease with increasing as predicted, however 
the detailed scaling is not clear. What is clear is that the radial propagation of the 
filament centre-of-mass increases with filament width and can outweigh the bi-normal 
displacement. 

4 . 3 . Comparison with the two-region model 

The two-region model, developed by Myra and Russell [371 EH], is an electrostatic 2D 
fluid model which includes both plasma resistivity and the effects of the magnetic 
geometry through a transformation of the cross-field differential operators. It has 
recently been extended to include the effects of warm ions jl3] . By including geometrical 
effects the two-region model is relevant for comparison with the simulations carried out 
here. The (cold ion) two region model has four regimes of filament dynamics which are 
characterized by the dimensionless parameters A and D where 



( 21 ) 


is a dimensionless resistivity and 



( 22 ) 


which represents the blob width. Each regime is characterized by a different 
closnre/dissipation mechanism for the diamagnetic (driving) current across the filament. 
The four regimes of dynamics, their current closure mechanism and their velocity 
scalings with Tg and S± are given in table For a derivation of the model and 
details see ref IS]. The resistive ballooning regime mentioned above is identical the 
the regime termed the inertially limited regime here. Both the Tg and scans trace 
specific paths in A, Q space. In calculating these paths the filament width has been 
converted from the definition used here (1% above background) to the Gaussian width 
of the filament, which is more appropriate for comparison with the two-region model. 
Geometric parameters are calculated directly from the flux-tube geometry. The least 
well constrained parameter required for this comparison is the parallel connection length, 
L||. In the two-region model L|| occurs when an integration of the curvature drive is 
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Table 1. Scaling of the normal velocity in the two-region model under differing regimes 
of current dissipation. 


Regime 

Closure Mechanism 

Velocity Scaling 

Resistive Ballooning (RB) 

Local polarization currents 

V oc 

Resistive X-point (RX) 

X-point polarization currents 

V cx 

Sheath Interchange (Cg) 

Sheath currents 

V oc 

Ideal Interchange (C/) 

Global polarization currents 

V oc 


carried out along the magnetic field line and balanced against either parallel or sheath 
currents. The curvature drive is assumed to be fixed at the midplane value such that 

2csPsi • Vln(n)dL t /23) 

Jo c)j_ 

In reality however the transformation of the filament cross-section due to flux expansion 
and magnetic shear can strongly enhance the curvature drive in regions where the 
perpendicular scale length of the filament becomes small. The transformation of the 
V operator along the flux tube is given by where the transformation matrix 

M. is given by [37] 

^ i^dBe^d) / {RuBo^u) {h — lu) RuRdBe,uBe,d/B 
0 (BuRuBe^u) / {BdRdBe,d) 

where R is the major radius, Bq is the poloidal magnetic field strength, B is the total 
magnetic field strength, / is the magnetic shear integrated along the magnetic field line 
and subscripts u, d refer to upstream (midplane) and downstream values respectively. 
The diagonal terms represent the decreased gradient in the normal direction and 
increased gradient in the bi-normal direction as a result of flux-expansion. The off 
diagonal term represents the increase in the gradient in the normal direction as a result 
of the shearing of the flux tube. The integration of the curvature drive in the two-region 
model can now be replaced by 

/ " 2csPs^ ■ V In (n) dL = 2csPs / " € ' Ad In (n) dL L^ff (25) 

Jo Jo Ox 





1 

1 


where 


= i,)M 


dL 


(26) 
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can be calculated directly from the flux-tube used in the simulations. Here gradient 
operator has been approximated by 

(27) 




l/5i 


/ 


When comparing with the two-region model, |^| L\\ can now be replaced with L^ff to 
provide a consistent comparison. Figure shows the scan paths using Le// in A, 
n space. Also shown is Le// compared to S, the distance along the flux tube from 
the midplane. The paths of the Tg and scan are both encompassed mainly in the 



Figure 14. Left; Le// calculated along the MAST SOL flux tube. Integration is 
carried out starting at the divertor. In the divertor region {S < —5ni) Le// begins to 
depart dramatically from L||. Right: Corrected paths in A, O space. 


RB regime, in agreement with the velocity and displacement scalings observed in the 
simulations. The suggests that, despite the complex nature of of the filament motion in 
3D, the velocity scaling normal to the flux surface may still be relatively well predicted 
by the two-region model so long as the effect of the flux-tube transformation on the 
curvature drive is accounted for. This comparison therefore supports the use of reduced 
models like the two-region model as a useful tool for comparison despite its neglect of 
3D physics, so long as suitable care is taken in determining Ln. 


5. Particle Transport 

Analysis of the centre of mass is an effective way of characterising the velocity scaling 
across parameter scans, but fails to capture the detailed transport of particles due to the 
evolution of the filament. A more robust analysis of particle transport can be carried 
out by calculating the particle flux crossing a set of concentric surfaces at each radial 
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location. This provides a radial profile of the instantaneous particle flux due to filament 
motion, , given by 

{x,t) = j Hf {x,t)v^ {x,t) dA{x) (28) 

Here x labels the normal coordinate, u/ is the filament density, is the E x B velocity 
in the normal direction and dA is the surface area of a surface located at a:. To iso¬ 
late the dynamics of the filament at the midplane each surface encompasses the entire 
toroidal domain but extends only 5mm either side of the midplane in the poloidal direc¬ 
tion; as such the surfaces resemble concentric toroidal bands centred on the midplane. 
By integrating these instantaneous particle fluxes in time it is possible to measure the 
total number of particles crossing each of the constructed toroidal surfaces during the 
evolution of the filament. This is shown in figure pT| for filament simulations across the 
temperature scan. 



Figure 15. Total number of particles crossing each radial location over the simulation 
time. Little variation is observed between different values of Tg, though a small fraction 
of particles reach a further radial extent at F = 20eV. The vertical dashed line 
indicates the starting position of the filament. 


The displacement of particles from their original position increases as the tempera¬ 
ture increases. This increase is weak however with the transport of the bulk of the 
particles remaining relatively constant between simulations. Indeed it is unlikely that 
variation in particle transport due to filament motion on the scale observed in figure 
may be observable above experimental uncertainty. This indicates that the role played 
by electron temperature in determining particle transport mediated by filaments may be 
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weak. Also notable is the fact that only a minority of particles are transported distances 
> hx • Particle transport into the outer region of the radial domain occurs in higher tem¬ 
perature filaments which are strongly turbulent. To explore whether this turbulence is 
driving this extra particle transport a second simulation at 20eV has been run with the 
Boltzmann source removed thereby removing drift-waves from the system. Comparing 
the two cases, as shown in figure [T^ it is clear that the particle transport towards the 
edge of the domain is a result of the drift-wave turbulence. Therefore, whilst drift-waves 
can act to disperse filaments [301 ED El they can also enhance particle transport beyond 
the levels expected in 2D theories. 



Figure 16. Comparison of the total number of particles crossing each radial position 
with (B) and without (NB) the Boltzmann source included. This indicates that the 
turbulence induced by the Boltzmann source can enhance particle spreading, though 
the enhancement is weak. 


In contrast to the scan in electron temperature, variation of the filament width, 
can result in a significant change to the particle transport within the filament. This is 
exhibited in figure [T7| where the 8cm filament shows the greatest spreading of particles. 
This is a combined effect of the increased size of the filament containing more particles 
and of the decrease in the deceleration of the filament with increasing S± (see figure [T^ . 
The magnitude of the different curves in figure [T^ is largely determined by the number 
of particles in the initial filament, which is clearly largest for the largest filament. The 
position of the peak of each curve however is determined by the motion of the filament, 
since if all filaments move at the same rate, then the surfaces that particles cross remain 
unchanged. Since the peak in the number of particles crossing radial surfaces is more 
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X (cm) 


Figure 17. Total number of particles crossing each radial location over the S± scan. 
A significant increase in the spreading of particles is observed as S± is increased. 

displaced for larger filaments this indicates that larger filaments are also more efficient 
at transporting particles into the far SOL. Coupling together the observations that large 
filaments both contain more particles and transport them more efficiently suggests that 
large filaments are likely to dominate particle transport into the far SOL. 

6. Discussion and conclusions 

This paper has sought to characterize the transport of particles at the outboard midplane 
in the SOL, using a geometry based on the MAST tokamak, due to filament motion 
through numerical simulation using the BOUT++ code. The presence of parallel density 
gradients, arising as a result of the magnetic geometry through ballooning of the filament 
and enhanced diffusion in the divertor region, drives the growth of a component of 
potential with even parity about the filament centre. This can cause a departure of the 
filament dynamics from standard 2D behaviour. These effects are most prominent in 
hotter filaments but begin to have an effect at low temperatures easily obtainable in 
SOL conditions. The presence of an even component of potential leads to an asymmetric 
potential dipole in the bi-normal direction causing rotation of the filament density before 
a break down into drift-wave turbulence occurs. Drift-waves are observed in all filaments 
across a temperature scan, however in cold filaments they appear as linear instabilities 
perturbing the filament front but not impacting the motion of the filament significantly. 
In hotter filaments the drift-waves are non-linearly driven and have a major impact on 
the filament motion. The large bi-normal displacement may suggest that inter-filament 
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interaction may occur, which cannot be captured by isolated hlament simulations and 
would require significantly more demanding simulations of turbulence. 

To reconcile the 3D filament simulations with 2D theories the scaling of the centre- 
of-mass velocity has been investigated across both a scan in temperature and filament 
width. The velocity of the filament in the direction normal to the flux surface show the 
characteristics of inertially limited filaments, or to use the terminology of the two region 
model, are in the inertially limited regime. The simulations were shown to be consistent 
with the two-region model only when the curvature drive in the divertor region was 
accounted for. The velocity of the filament in the bi-normal direction showed several 
features of a scaling based on the Boltzmann response, however the complex dynamics 
of the filament in this direction stretch the applicability of a velocity scaling argument. 
It is interesting that the normal and bi-normal velocities appear to show a degree of 
de-coupling, and that 2D models may still be able to capture the motion in the normal 
direction relatively well. This is an important point because the most important aspect 
of particle transport is that which is occurring in the normal direction. If reduced 2D 
models are capable of capturing the salient features of transport in this direction, despite 
the neglect of 3D effects which strongly perturb the dynamics of the filament, then 
they may be a valuable and efficient tool in determining radial particle transport due to 
filament motion. It is notable that the displacement of the centre of mass of the filaments 
simulated never exceeds a filament width. In experiment on MAST filaments are 

routinely observed to propagate further than their width. This apparent disagreement 
between simulations and experiment may be rooted in the fact that the model used is 
a physically reduced system alongside the difficulty in interpreting experimental data. 
In terms of the model, it is possible that the addition of extra physics into the system 
will provide more coherent filaments. Finite ion temperature has been shown to provide 
such an increase in coherence for smaller blobs [M] whilst electromagnetic effects have 
recently been shown to increase the coherence of high /5 filaments m- Neither of these 
effects are included here and may be pursued in the future, though given that filaments 
in MAST are generally larger than the gyroscale and are relatively low pressure it is 
unclear how important these two effects may be. Furthermore whilst the centre of 
mass of the filament in simulations propagates a limited distance the density peak of 
the filament propagates much further. This presents some ambiguity when comparing 
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the propagation of filaments in experiment with simnlations. This work would greatly 
benefit from a careful validation exercise between simulation and experiment. Such an 
exercise is already planned for MAST in the near future. 

It is worth noting that the filament amplitude has not been varied in this study, however 
in reality will provide another mechanism by which the filament velocity may change. 
Angus et. al [53 have recently demonstrated that for filaments with amplitudes of unitary 
magnitude, the Boussinesq approximation can break energy conservation which can 
impact the dynamics of the filament. For large amplitude filaments this effect can 
be dramatic, however for unit amplitude filaments the deviation is less significant. 
Furthermore the scalings with respect to temperature and filament width are not 
likely to change significantly as a result of the Boussinesq approximation. It may be 
possible however that the cause of the deviation of these simulations from experimental 
observations on MAST is the amplitude of the simulated filaments. In experiment large 
amplitude filaments exist, and if it is these filaments that propagate into the far SOL, 
then they must be modelled. If this is to be achieved then the full system of equations 
must be considered without the Boussinesq approximation, perhaps in the manner of 
Angus et.al |55] - 

A deeper analysis of particle transport at the outboard midplane has been achieved 
by calculating the particle flux induced by the motion of the filament through a 
series of concentric toroidal bands centred on the midplane. This has shown that the 
presence of drift-wave turbulence in hotter filaments provides a slight enhancement to 
particle transport by comparison between two simulations, one with and one without 
the Boltzmann source included (which drives the turbulence). This enhancement is 
small however and, across the scan in the particle transport due to the filament 
motion remains broadly similar. By contrast the filament width plays a major role in 
determining the radial transport of particles, with particles in larger filaments being 
transferred much further in the normal direction. This suggests that a dominant 
parameter determining the filamentary contribution to radial particle transport is the 
distribution of filament widths as opposed to background plasma temperature. 
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